Nonlinear Volatility of River Flux Fluctuations 
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We study the spectral properties of the magnitudes of river flux increments, the volatility. The 
volatility series exhibits (i) strong seasonal periodicity and (ii) strongly power-law correlations for 
time scales less than one year. We test the nonlinear properties of the river flux increment series by 
randomizing its Fourier phases and flnd that the surrogate volatility series (i) has almost no seasonal 
periodicity and (ii) is weakly correlated for time scales less than one year. We quantify the degree 
of nonlinearity by measuring (i) the amplitude of the power spectrum at the seasonal peak and (ii) 
the correlation power-law exponent of the volatility series. 
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Climate is strongly forced by the periodic variations 
of Earth with respect to state of the solar system. The 
seasonal variations in the solar radiation cause to pe- 
riodic changes in temperature and precipitation which 
eventually lead to seasonal periodicity of river flow. In 
spite of this well defined seasonal change, river flow ex- 
hibits highly unpredictable complex behavior; floods and 
droughts are usually unexpected and cause severe dam- 
age in life, housing, and agriculture products. Hence, 
river flow is likely to have indirect nonlinear response to 
the seasonal changes in solar radiation. 

Many components of the water budget of a catchment 
are coupled in a nonlinear fashion. The key for all in- 
teractions between atmospheric processes like precipita- 
tion, temperature, humidity (or extraterrestrical inputs 
like sun radiance) and surface runoff is the soil. The 
dynamic state of this key variable is highly nonlinear. 

By means of the proposed methods in this paper it 
will be possible to characterize quantitativly the degree 
of nonlinearity of the involved processes by investigating 
the outputs of the catchment (the resulting flux time se- 
ries) only. This nonlinearity test would be very helpful 
both for the design of time series models and statistical 
prediction algorithms. 

There are several statistical features found in the ear- 
lier studies of river flow fluctuations. E.g., river flow fluc- 
tuations have broad probability distribution [1,2]. More- 
over, river flow fluctuations have unique temporal orga- 
nization; they are long-range power-law correlated and 
possess scale invariant structure [3,4]. These river flow 
power-law correlations are usually characterized by scal- 
ing exponent [5,6] as was originally defined by Hurst in 
his seminal work [3] regarding the Nile river floodings. 
However, such scaling laws only quantify the linear prop- 
erties (two-point correlations) of a time series. Here we 
study other nonlinear statistical aspects of river flow fluc- 
tuations. 

A nonlinearity of a stationary time series may be de- 



fined with respect to its Fourier phases [7,8]. Series that 
its statistical properties are independent of the Fourier 
phases may be defined as linear otherwise the series 
may be defined as nonlinear. Autoregression processes 
and fractional Brownian motion are examples for lin- 
ear processes while multifractal processes are examples 
for nonlinear processes. Recently it has been shown 
that volatility correlations of long-range power-law cor- 
related time series reflects the degree of nonlinearity of 
a time series [8]. Given a time series, Xi, the volatility 
series is defined as the magnitudes of the series incre- 
ments, [Axi] = \xi+i — Xi\. It was found that long-range 
correlated linear series has uncorrelated volatility series 
while long-range correlated nonlinear series has corre- 
lated volatility series; see [8] for details. Power-law cor- 
relations in the volatility series indicate that the mag- 
nitudes \/S.Xi\ are clustered into self-similar patches of 
small and big magnitudes — a big magnitude increment 
is likely to precede a big magnitude increment and vice 
versa. When the volatility series, [Axi], is uncorrelated 
the increment series is homogeneous. Volatility correla- 
tions were found, for example, in econometric time series 
[9], heart interbeat interval series [8,10], and human gait 
dynamics [11]. 

Here we study the volatility properties of river flow 
fluctuations. We first extend the notion of volatility 
to periodic time series by measuring the periodicity of 
the volatility series. We find that after randomizing the 
Fourier phases of the river flow increment series, the pe- 
riodicity of the volatility series is almost diminished in- 
dicating that "periodic volatility" is a result of nonlin- 
earity. We also find long-range volatility correlations for 
time scales below one year. Our results suggest that clus- 
ters of magnitudes of river flow increments appear in two 
ways: periodic clustering and long-range self-similar clus- 
tering. 

We analyze the daily river flux time series of 30 world 
rivers scattered around the globe. The mean flux of these 



1 



rivers ranges from ^ 0.6m^/s to ^ 2 x lO^m^/s and thus 
covers more than 5 orders of magnitudes. The length of 
the these time series ranges from 26 years to 171 years 
with average length of 81.3 years. In Fig. 1 we present a 
typical example of 4 years (1986-1990) of river flow data 
of the Maas river in Europe. It is evident that fluctu- 
ations around large river flow are also large while the 
fluctuations around small river flow are small. 

To study the nonlinear properties of the river flow time 
series we apply a surrogate data test to the river flow 
increment series. We use a surrogate data test that pre- 
serves both the power spectrum and the probability dis- 
tribution of the river flow increment series [7]. On the 
other hand, the Fourier phases of the surrogate series 
are random. Thus, the surrogate data test linearizes the 
series under consideration. Since the histogram of the 
surrogate data is identical to the histogram of the origi- 
nal increment series one can be sure that the probability 
distribution is not the source of the nonlinearity of the 
data. Fig. 2 shows the river flow increment scries and 
its power spectrum before and after the surrogate data 
test. Although the river flow increment series exhibits 
irregular behavior, its power spectrum shows a very pro- 
nounced seasonal peak with few harmonics. As expected, 
the surrogate series shows a similar pattern with almost 
identical power spectrum. 

Next we compare the power spectrum of the volatil- 
ity series obtained from the original increment river flow 
series and from the surrogate scries (Fig. 3). The 
power spectrum of the original volatility series shows a 
pronounced seasonal peak while the power spectrum of 
the surrogate volatility series has no seasonal periodic- 
ity. The seasonal periodicity of the original volatility 
series may be associated with the increased fluctuation 
for large river flux (see Fig. 1). The absence of seasonal 
periodicity for the surrogate volatility series is somehow 
counter intuitive since the surrogate series itself is as pe- 
riodic as the original river flow increment series while a 
simple inversion operation of the negative values of Axi 
to obtain |Aa:,| diminishes this periodicity. The absence 
of the seasonal periodicity from the surrogate volatility 
series indicates that periodicity in the magnitude series 
is a result of nonlinearity. We suggest that the ampli- 
tude of the seasonal peak of the original volatility series 
compare to the seasonal peak of the surrogate volatility 
series would quantify the degree of nonlinearity. 

We use the power spectra of the original and surro- 
gate volatility series to analyze the correlations proper- 
ties of these series. If a series Xi is long-range correlated 
than its autocorrelation function decays as a power law 
C{1) = J2i=i ^i+i^i ~ where A'' is the series to- 
tal length, I is the lag, and 7 is the correlation exponent 
(0 < 7 < 1). In this case also the power spectrum fol- 
lows scaling law [12,13] S{f ) ~ l/f>^ where 7 = 1 - In 
Fig. 4 we show the power spectra of the original and sur- 
rogate volatility series for frequencies larger than lyr^^ . 
There is a notable difference between the power spectrum 



of the original volatility and the power spectrum of surro- 
gate volatility; while the power spectrum of the surrogate 
volatility series is almost flat, the power spectrum of the 
original volatility decays as a power law with an expo- 
nent of /3 « 0.66. Thus (i) the original volatility series is 
power-law correlated and (ii) its correlations are a non- 
linear measure since they significantly reduced after the 
surrogate data test. The interpretation of these correla- 
tions is that there are clusters of big magnitudes |AFj| 
that are statistically followed by patches of big magni- 
tudes. These clusters are in addition to the periodic 
clustering (shown in Fig. 3). We note that the power 
spectrum is not the preferred method for scaling analy- 
sis; we repeated the scaling analysis with more advanced 
method, the detrended fluctuation analysis [14] , and find 
less noisy but similar results [15]. 

We summarize the periodic volatility and the volatil- 
ity correlations results for the 30 rivers under consider- 
ation in Fig. 5. To systematically compare the seasonal 
periodicity of the different rivers we first normalize the 
volatility series by subtracting its mean and dividing by 
its standard deviation; thus, the area under the power 
spectrum of the different volatility series should be the 
same. The seasonal peak of the volatility series exists for 
all 30 rivers and is significantly higher then the seasonal 
peak of the surrogate volatility (Fig. 5 upper panel). The 
scaling exponent /3 of the original volatility series (Fig. 5 
lower panel) indicates correlations; in most of the cases 
(27/30=90%) the exponent of the original volatility se- 
ries lies above 1 standard deviation of the exponent of 
the surrogate volatility scries. The average ± 1 standard 
deviation of the scaling exponent of original volatility se- 
ries is /? = 0.49 ±0.11 and is significantly higher than 
the average ± 1 standard deviation of the scaling expo- 
nent of the surrogate volatility scries (3 = 0.18 ± 0.13; 
the p value of the Student's t— test is less than 10^^. For 
time scales larger than 1 year the volatility series is only 
weakly correlated with average exponent f3 = 0.27 ±0.26. 

Thus, we find two measures of nonlinearity related to 
the river flow data, periodic volatility and long-range cor- 
related volatility. These two measures are related to the 
clustering of the magnitudes of river flow fluctuations, 
periodic and long-range correlated clustering. 

To study in more details the possible source for such a 
seasonal periodicity of the volatility series we propose a 
simple scheme to generate series with some similar char- 
acteristic properties as in the river flow data. To mimic 
the enhanced fluctuations for large river flow we assume 
that, 



Xi = {1+ Ar]i)si 
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where ry, is a Gaussian random variable with zero mean 
and unit standard deviation, A is the noise level, and Sj 
is an asymmetric periodic function. 



Sj+nT — 



_ / 1 + cos(27r/i) for < j < |T 
l-cos(47r/i) for |T < j < T, 



(2) 



2 



where T is the time period T = 365 in arbitrary time 
units, j is an integer < j < T, f = 0.75/r, and n 
is an integer. Xi decreases for 2/3 of the time period 
T and increases for 1/3 of this time period. When the 
noise level A increases the nonlinear term ArjiSi also in- 
creases. We generate Xi series with different noise level 
and then calculate the power spectrum of the normal- 
ized volatility series |A.7:,;| of the original and surrogate 
Aa;^ scries (Fig. 6). We find that when the noise level is 
relatively small the seasonal peak is present in both the 
original and surrogate volatility series. When the noise 
level increases the periodicity of the surrogate volatility 
series is diminished. Thus, the larger is the difference 
between the peak of original volatility series and peak 
of the surrogate volatility series, the larger is the non- 
linearity of Axi. This simple scheme indicates that the 
surrogate data test not always deminishes the seasonal 
periodicity of the volatility series but rather eliminates 
the nonlinear part of the process that is proportional to 
the noise level (Eq. 2). We also analyzed time series 
generated by realistic hydrological model (ASGi model 
[Kontinuierlicher Abfluss und Stofftransport- Integrierte 
Modellierung unter Nutzung von Geoinformationssyste- 
mcn] for Bavaria, Germany [16]) for three rivers: Naab, 
Regniz and Vils. Both, the seasonal periodicity of the 
volatility series as well as its correlations are reproduced 
by the model and disappear after phase randomization, 
as was observed in the real data. 

In summary we analyse the periodic volatilities and the 
time correlations of river flow data for 30 world rivers. 
We find that the volatility series are strongly correlated 
with a power law behaviour for time scales less than 1 
year. The periodic volatility and the long-range corre- 
lated volatility are found to disappear when randomiz- 
ing the phases. This indicate that these features of the 
volatility time series are due to nonlinear dynamical pro- 
cesses. We suggest that such nonlinear features may re- 
sult from an interaction between noise and the seasonal 
trends. 

Preliminary analysis of other climate records, such as 
daily temperature and daily precipitation records, shows 
the existence of periodic and long-range volatility with 
similar properties as for the river flow data. Thus, the 
results presented here may be generic for other climate 
records. 
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FIG. 1. Typical river fiow time series of Maas river (Eu- 
rope, 1986-1990). The record shows a periodic pattern with ir- 
regular fluctuations. Note the large fluctuations around large 
river fiow compare to small fluctuations around small river 
flow. 
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FIG. 2. River flux increment series of Maas river (left pan- 
els) and their corresponding power spectra (right panels) be- 
fore (upper panels) and after (lower panels) the surrogate test 
for nonlinearity. The series length is 80 years where just the 
last 4 years of the record are shown in the left panels. The 
original river flow increment series and the surrogate incre- 
ment series have identical probability distribution and very 
similar power spectrum. 
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FIG. 3. Same as Fig. 2 but for the river flow volatility 
series, |AFi|. Here, the original volatility series shows a pro- 
nounced seasonal peak while the surrogate volatility series 

doesn't show such a peak indicating that the periodicity in 
the volatility series is a result of nonlinearity. 
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FIG. 4. Log-log plot of the power spectra shown in Fig. 3. 
The solid lines are the best fit of S{f) ~ 1//'^ for frequen- 
cies 1.05j/r~^ < / < 52yr~^; we use logarithmic binning for 
the exponent calculation. The original volatility series (left 
panel) decays as a power law i/j''-" *'*' indicating long-range 
correlations of the volatility series. The power spectrum of 
the linearized surrogate volatility series has a flater spectrum 
indicating much less correlated behavior. Thus, correlations 
in the volatility are an additional measure for nonlinearity of 
the river flow increment time series. The vertical dashed lines 
indicate 1 week periodicity. 
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FIG. 5. A summary of the results obtained for the 30 world 
rivers. For each river flow increment series (•) we generated 10 

surrogate series (o) and calculated the amplitude of the sea- 
sonal peak of the volatility series (upper panel) and the scal- 
ing exponent /? for frequencies 1.05yr~^ < / < 52yr~^ (lower 
panel); the average and 1 standard deviation are shown. In 
order to systematically compare the results of the different 
rivers we subtract from the volatility series its mean and nor- 
malized by its standard deviation. The seasonal peak of the 
volatility series is significantly higher compare to the seasonal 
frequency of the surrogate volatility series (upper panel). The 
scaling exponent /3 shown in the lower panel is systematically 
higher for the original volatility series. For 28 rivers the orig- 
inal volatility exponent is larger than its surrogate exponent 
where 27 of these 28 exponents lie well above the error bars. 
The error bars on the right hand side are the group average ± 
1 standard deviation for the original and surrogate volatility 
scaling exponent. 
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FIG. 6. The power spectrum of the normahzed volatUity 

series, |Aa::i|, of an artificial time series a;^ = (1 + Arii)si for 
different noise level A; sec Eqs. (1),(2). The upper pan- 
els shows the power spectra of the original volatility series 
and the lower panels show the power spectra of the surrogate 
volatility series. When the noise level increases (from the left 
panels to the right panels) the seasonal peak of the surrogate 
volatility series reduces. The harmonics of the power spectra 
are partly caused by the asymmetric Xi and partly because of 
the absolute value operation for the volatility series. 
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